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Measuring the small primordial nonGaussianity (PNG) predicted by cosmic inflation theories may 
help diagnose them. The detectability of PNG by its imprint on the 21cm power spectrum from 
the epoch of reionization is reassessed here in terms of /jvl, the local nonlinearity parameter. We 
find that an optimum, multi-frequency observation by SKA can achieve A/jvl 3 (comparable to 
recent Planck CMB limits), while a cosmic- variance- limited array of this size like Omniscope can 
even detect AfNL ~ 0.2. This substantially revises the methods and results of previous work. 
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Introduction. — The theory of cosmic inflation [J Q, 
advanced to solve the cosmological horizon and flatness 
problems, also explains the initial fluctuations which 
later gave rise to galaxies and large-scale structure in the 
universe. While inflation generically predicts initial mat- 
ter density fluctuations with an approximately Gaussian 
random distribution, the small deviations from Gaussian- 
ity that characterize different inflation models has been 
suggested to provide an observational probe to test and 
distinguish the models. While purely Gaussian initial 
density fluctuations are fully described by their power 
spectrum, primordial non-Gaussianity (PNG) requires 
higher-order statistics to characterize it, the lowest-order 
being the 3-point correlation function, or its Fourier 
transform - the bispectrum - which is zero for the Gaus- 
sian case. Henceforth, we will describe the level of PNG 
predicted by different inflation models in terms of this 
bispectrum as parametrized by the dimensionless "non- 
linearity parameter" f]^'^^, specialized to the case of the 
so-called "local" template. [For further details, see [sl 
and refs. therein.] 

The standard simplest model - slow roll, single-fleld in- 
flation - predicts an extremely small level of PNG, given 
by f]^}^^ = (5/12)(l — Ug), where rig is the spectral index 
of the primordial power spectrum 0-0] • Recent cosmic 
microwave background (CMB) temperature anisotropy 
measurements find ris « 0.96 so /j^^'^' « 0.016. 

Other more general models (e.g. multi-fleld inflation) pre- 
dict much larger values of /]vF'- 

Observational cosmology has made important progress 
in constraining PNG thus far. Recent measurement of 
the CMB anisotropy bispectrum by Planck has placed the 
most stringent constraint so far, f]^}^^ = 2.7 ± 5.8 [11). 
Since, even in the ideal noise-free limit, CMB tempera- 
ture (temperature-|-polarization) measurements can only 
reduce the error to A/]^f ^ « 3.5 (1.6) [B, Elj, there is 
great interest in flnding other methods to measure PNG; 



if future observations still do not detect PNG, an error 
budget A/]^£^' < 1 will be necessary to rule out non- 
standard inflationary models conclusively. [Henceforth, 
we focus on this "local" template and remove the label 
"local".] 

PNG affects the clustering of the early star-forming 
galactic halos responsible for creating a network of ion- 
ized patches in the surrounding intergalactic medium 
(IGM) during the epoch of reionization (FOR), which 
leaves a PNG imprint on the tomographic mapping of 
neutral hydrogen in the IGM using its redshifted 21cm ra- 
diation. We shall here investigate in detail the prospects 
for constraining PNG with radio interferometric 21cm 
measurements. Our method and results differ signifi- 
cantly from previous attempts in the literature 



as follows: (1) We apply the ionized density bias derived 
by Ref. Q to model the effect of PNG by the excursion- 
set model of reionization (ESMR); (2) We show a phe- 
nomenological model that can constrain PNG just as ac- 
curately, independent of reionization details; (3) We show 
that a single-frequency (i.e. single epoch) measurement 
can be tuned to the optimum frequency for constrain- 
ing PNG; (4) We show how combining multi-frequency 
(multi-epoch) measurements further reduces the forecast 
errors. 

PNG Signature in the 21cm power spectrum. — The 3-D 
power spectrum of 21cm brightness temperature fluctu- 
ations (hereafter, "21cm power spectrum") in observer's 
redshift space can be expressed to linear order in neutral 
and total hydrogen density fluctuations, Sp^^ and Sp^ , re- 
spectively, as the sum of powers of = k • n/[k[ (cosine 
of angle between line-of-sight (LOS) n and wave vector 
k of a given Fourier mode) [l^ Il3li Pat(M,z) = 

Ps,^, .5pHi ik,z) + 2 Ps^^^ ,s,^ {k, z) fil 
ik,z)^l , where Sfb{z) = (23.88 mK) 
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FIG. 1. The spherically-averaged 21cm power spectrum at 
z — 10.10, using ESMR with the values of (/jvl,Cesmr) 
marked in the legend. 
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QmIi^ TcT' ^'^"^ ^Rii^) is the global neutral fraction. 
Pa.b is the power spectrum between fields a and b. 
Here, we focus on the limit where spin temperature 
Ts TcMB, valid soon after reionization begins. As 
such, we can neglect the dependence on spin tempera- 
ture, but our discussion can be readily generalized to 
finite Ts. We also focus on the 21cm signal on large 
scales k < 0.15 Mpc~^, so that linearity conditions are 
met (see Ref. 18] for a summary of these conditions). 
When the typical size of ionized regions is much smaller 
than the scale of interest, nonlinear effects of reionization 
patchiness on the 21cm power spectrum [3| can be 
neglected. If we define neutral and ionized density 
biases, bp^^ and fopnn, according to ba{k) = Sa(k-)/Sp(k), 
i.e. ratio of density fluctuation in field a to that of total 
matter density in Fourier space, then the 21cm power 
spectrum can be rewritten as 

. 2 2 

FAT(k,z) = (5r^5|i [6pHi(fc,2) + Atk] Pssik,z), (1) 

where Pssik,z) is the total matter density power spec- 
trum. Here, we assume the baryon distribution traces 
the cold dark matter on large scales, so Sp^ = 5 p. 

The ionized density bias is the fundamental quantity 
derived from reionization models, related to the neutral 
density bias by 



bpui = (1 - a^HiifepHiJAffl 



(2) 



where xhii = 1 — xm- We model reionization with PNG, 
as follows, based on the results of Ref. Q. 

(i) ESMR: The basic postulate of ESMR[i3| is that 
the local ionized fraction within a spherical volume with 
radius R is proportional to the local collapsed fraction 
of mass in luminous sources above some mass threshold 

Mmin, i.e. Xmi{M^in, R, z) = CeSMR /coll(-Mmin, i?, z), 

where Cesmr parametrizes the efficiency of this mass in 
releasing ionizing photons into the IGM. For simplicity, 
we assume atomic-cooling halos (ACHs), i.e. halos with 



virial temperature Tvir ^ 10^ K, are the only sources of 
ionizing radiation. 

Our methodology for ESMR with PNG is as follows. 
The collapsed fraction of ACHs in Ref. [i^l (see also 



21|). calculated in the non-Mar kovian extension to the 
excursion set formalism 22-23|, for a given PNG pa- 
rameter, is applied to the ESMR formalism to calculate 
the ionized density bias, for a given efficiency, analyt- 
ically, as described in detail in Ref. Q. Henceforth, 
since the functions xmi{,z) and bp^^^{k,z) are set by 
two parameters (/atl, Cesmr), given a fiducial cosmol- 
ogy, so is the 21cm power spectrum PAT(k, 2) at any 
z. As Figure [T] illustrates, while Cesmr changes the am- 
plitude of the 21cm power spectrum, /jv^ changes the 
shape at small k significantly. The (nonzero) minimum 
of the curve for /jvl < is at wavenumber k-^^ where 
the non-Gaussian 6pjjj(fci,z) ss —1/3, since Eq. ([T]) gives 

(ii) Phenomenological ("Pheno"-) model: Just as PNG 
exhibits a scale-dependent effect on halo bias 23, 21. 25- 
1301] . so we also find, in Ref. 01, a scale-dependent 
non-Gaussian correction to the ionized density bias, 
A5j,Hn(fc,z) = 6pH„(A:,z) - b%^^{z), where 6^^^^ is the 
Gaussian ionized density bias and scale-independent. 
[According to Ref. Q, there is also a scale-independent 
non-Gaussian correction, A6pHjj(z). However, Afop^jj <C 
^/?Hii ^'^'^ f'^L < 10, so we neglect it here, similar to the 
neglect of a scale-independent non-Gaussian correction to 
the halo bias when constraining PNG with galaxy surveys 

For the local template, we derived from the ESMR a 



relation between A6pjjjj and b'^^ 
scales, 



in Ref. [3]. On large 



Ae(^.=)-v„.[c,w-i]J2|i|^, (3) 

where Sc « 1.686 is the critical density in the spherical 
collapse model (in an Einstein-de Sitter universe), D{z) 
is the linear growth factor normalized to unity at z = 0, 
5(0) = (1 -I- Zi)~^D~^{zi) « 0.76 in our fiducial cosmol- 
ogy, where Zi corresponds to the initial epoch, i.e. limit 
of large redshift, and T{k) is the matter transfer function 
normalized to unity on large scales. 

This relation was further tested and confirmed by 
numerical solution of the linear perturbation theory of 
reionization (LPTR) which includes radiative transfer 
[3^ . [For further details, see Ref. [1].] Henceforth, 
Eq. (j3|) is assumed to be generic, regardless of reioniza- 
tion details. In what we call the "Pheno-model" , the 
21cm power spectrum is set by three parameters - /nl, 
xhi and b'^^^^, at a given redshift. The latter two param- 
eters embrace our ignorance of reionization. 

Observability of 21cm interferometric arrays using Fisher 
matrix formalism. — Radio interferometric arrays mea- 
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TABLE I. Specifications for 21cm interferometers. We as- 
sume observation time to = 4000 hours for each redshift bin 
of bandwidth B = 6 MHz. 



TABLE II. Forecast l-a errors at xmi ~ 0.25(2 = 11.24) and 
xnii ~ 0.50 {z — 10.10), respectively. "F.V." means "fiducial 
values" . 



Experiment Nin Lmin (m) 77 ^e(2: = 6/8/12)[m^] 



MWA^ 


50 


12.5 


1 


9/14/18 


LOFAR 


32 


100 


0.8 


397/656/1369 


SKA 


1400 


10 


0.8 


30/50/104 


Omniscope 


10^ 


1 


1 


1/1/1 



»[sr] ^ 

2(AV^,X 
2tt 



^ Sky rotation adds an additional factor of two to tlie actual 

observed patches. 
''128 antennae total [3^. 

Assume LOFAR can simultaneously observe two patches on the 

sky. 



sure the 21cin signal from coordinates = OxS-x ^(^y^y + 
Az^n, where {6x,0y) mark the angular location on the 
sky, and Ai> is the frequency difference from the cen- 
tral redshift z^, of a redshift bin. It is related to the 
3D Cartesian coordinates r (with origin at the bin cen- 
ter) by ©j_ = r±/rfyi(2:*), and Av — r\\/y{zt), where 
dA{z) is the comoving angular diameter distance, y{z) = 
\2i{l + zf/H{z), A21 = \{z)/{l + z) w 0.21m, and H{z) 
is the Hubble parameter at z. The Fourier dual of is 
defined as u = u^fix + UyCy -f M||n (u|| has units of time), 
which is related to k (Fourier dual to r) by Uj_ — (i^kj_ 
and W|| = y The power spectrum in u-space is related 
to that in the k-space by Pat(u, z) — PAT(k, z)/{d\y). 

For an interferometric array, a baseline L corresponds 
to — 27rL/A. (Note: our convention is differ- 
ent from the observer's convention uj_ = L/X.) Let 
n(Lu^)d^L denote the number of redundant baselines 
Lu^ corresponding to uj^, i.e. autocorrelation of ar- 
ray density. Then, the noise power spectrum in u- 



space 

[3 



1,111 is P^M = (^) / M(LuJ], where 

« (280 K) [(1 -I- z)/7A]^'^ is the system temperature 

cx (for A less than antenna size) is effective 

collecting area, and tg is total observation time. 

We adopt the following configuration of interferometric 
arrays. We assume antennae are concentrated within a 
nucleus of radius i?o with area coverage fraction close to 
100%, with coverage density dropping like in a core 
extending from Rq to iij^. We neglect the dilute antennae 
distribution in the outskirts, R > Ri^. Given central 
array density po, the configuration can be specified by 
two convenient parameters: Mn, the number of antennae 
within i?in, and 77, the fraction of these antennae that 
are in the nucleus . The relations are Rq = y^rjNin/ pqtt, 
i?i„ = i?oexp[(l-?7)/(277)] H. 

Given a parameter space {pa}, the Fisher matrix 
for 21cm power spectrum measurements is Fat — 
(^^) (^^) / [SPATin)f. The l-a fore- 
cast error of the parameter pa is given by Apa = 
•\/ (r~-'-)aa- The power spectrum measurement error in 
a pixel at u is SPat{u) = [Pat{u) Pn{u±)] / 



ESMR 



Pheno-model 



a;Hii 






/nl 


Cesmr 


/iVi 




PHII 






[F.V. 





50.0 





0.75 


6.19 J 




MWA 




13000 


1500 


14000 


300 


8800 


0.25 


LOFAR 




1200 


130 


1200 


1.1 


29 




SKA 




16 


1.8 


16 


0.028 


0.79 




Omniscope 




0.38 


0.040 


0.38 


0.00044 


0.012 






[F.V. 





50.0 





0.50 


5.43 J 




MWA 




700 


63 


750 


17 


220 


0.50 


LOFAR 




100 


8.1 


96 


0.16 


2.0 




SKA 




19 


1.5 


18 


0.030 


0.37 




Omniscope 




1.8 


0.15 


1.8 


0.0023 


0.027 



where Nc — uj_du±du\\ilB / {27r)'^ is the number of in- 
dependent modes in that pixel. We adopt logarithmic 
pixelization, du±/u±_ = (iu||/u|| = 10%. Here, D, is solid 
angle spanning the field of view, B is bandwidth of the 
redshift bin. 

We assume experimental specifications in Tab le m for 



MWA[36|, LOFAR137|, SKA[38l, and Omniscope [391, re- 
spectively. We assume residual f oreg rounds can be ne- 
glected for A; y > fc||_inin = 2TT/{yB) |34| . The minimum fc_L 
is set by the minimum baseline, fc_L,inin = 27rI/i„m/(AdA)- 
We account for modes up to fcmax = 0.15Mpc~^. Our 
fiducial cosmology is as follows: i}\ = 0.72, Hm ~ 0.28, 
= 0.046, Ho = lOO/ikms-^Mpc"^ {h = 0.7), erg = 
0.82, ris = 0.96, consistent with WMAP7 results [40|, 
with linear matter power spectrum of Ref. [4lj. ESMR 
fiducial values are /nl = 0, Cesmr = 50.0, corresponding 
to electron scattering optical depth Tcs = 0.08. To facil- 
itate direct comparison with ESMR, the fiducial model 
of 21cm power spectrum for the Pheno-model will be the 
same as the ESMR fiducial model, which sets the fiducial 
values of xhi and b'j^^^^ for a given redshift. 

Results. — 

(i) Single epoch constraints: In Table [Hi we list 
forecast errors, A/^vl, marginalized over Cesmr in the 
ESMR, and over xhi and b^^^^ in the Pheno-model. Top 
and bottom sets of forecasts use information from a sin- 
gle redshift bin centered at 2; = 11.24 (xhii = 0.25) and 
z = 10.10 (xhii = 0.50), respectively. For the same red- 
shift and experiment, the values of A/nl match very 
well between these two models. This demonstrates that 
the Pheno-model, which makes no assumptions about the 
connection between xm and ^p^^j , can constrain f^L for 
a single epoch measurement as accurately as a reioniza- 
tion model which links the evolutions of a;Hi and 6^^^^^ . 

For the same experiment, as Tables HI] shows, values of 
AJnl differ significantly between xhii = 0.25 and 0.50, 
implying a strong redshift dependence for A/tv^. For 
Omniscope, in particular, Afj^L shrinks by a factor of 
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FIG. 2. 1-cr error A/jvl from a series of single redshift bins, 
each with 6 MHz bandwidth. 



TABLE III. The 1-a error A/jvl, marginalized over ^esmr, 
from information that combines multiple redshift bins, each 
with 6 MHz bandwidth, ["l-band", "3-band", "5-band", "7- 
band"] means the result from [only the z = 11.24 bin {xmi = 
0.25), three redshift bins {xmi ~ 0.17—0.37), five redshift bins 
{xuii « 0.11 - 0.52), seven redshift bins (xmi ^ 0.06 - 0.70)], 
respectively. 



^max 


Experiment 


l-band 


3-band 


5-band 


7-band 




MWA 


13000 


1800 


520 


200 


0.15 Mpc"^ 


LOFAR 


1215 


91 


39 


26 


(default) 


SKA 


16 


5.0 


3.5 


2.8 




Omniscope 


0.38 


0.23 


0.18 


0.16 




MWA 


19000 


2100 


700 


280 


0.10 Mpc"^ 


LOFAR 


2600 


130 


62 


43 




SKA 


27 


8.1 


5.8 


4.8 




Omniscope 


0.84 


0.52 


0.42 


0.36 



~ 5 from xmi = 0.50 to 0.25. To investigate this in detail, 
we use the ESMR to plot A/jvl vs z in Figure [5] For 
MWA and LOFAR (both noise-dominated experiments) , 
the constraint is tighter at lower redshifts where the noise 
is smaller. However, for SKA and Omniscope (cosmic- 
variance-dominated experiments), there appears to be a 
"sweet spot," where A/atl is minimized, at xhii ~ 0.25 — 
0.30. 

This sweet spot for SKA and Omniscope 
can be explained using the Pheno-model, as 
follows. For the sake of argument, assume 
the ideal noise- free limit. The derivative 

\dPAT/dfNL\ oc (xmi/xm) [b%^, ~ l] \b%^ 



|(1 - 6g^J(i -f 6^^J|, using Eq. (0), while the power 
spectrum error APat fx Pat oc [6^^^ -I- /i^ ^ 
^ ■ ^ ' ■ ^ when averaged over /ik . 



"phi 



I ■'^hen averaged over /ik. Since 

Aa/a ^ 1 for a = zhi and b'^^^^ in Table [TT] for this 
limit, we can fix xm and b^^^^, so the Fisher matrix 

F = (A/A,i)-2 = [(dPAT/d/jvi) /APat]' is peaked 
when bf^^ « -(6 + \/2T)/15 « -0.71. Eq. ([2]) gives 
iffli (1 - ^pHi)/(^pHii ~ ^PMi)- ^^^^'^ Gaussian bias 
^Pmi ~ ^ (^"^^ Fig. 5 of Ref.[3|), the sweet spot is at 
Xhii ~ 0.25 in the ideal noise-free limit. In reality, 
finite noise in a cosmic-variance dominated experiment 
is larger at higher redshifts, so the sweet spot will occur 
at slightly lower redshift, depending on the amplitude of 
noise. 

(ii) Multi-epoch constraints: While a futuristic sin- 
gle epoch measurement can achieve a remarkable accu- 
racy of A/jvL = 0.38 (16) for Omniscope (SKA), adding 
tomographic information can further improve the accu- 
racy. The Pheno-model alone cannot be used to com- 
bine multi-frequency measurements because it does not 
specify the redshift evolutions of xhi and 6pjjjj. On the 
other hand, a model such as the ESMR can be used to 
combine multi-frequency measurements because it fixes 
the reionization history (and therefore Pat) for a given 
{/nl, Cesmr). We show multi-epoch constraints from the 



ESMR in Table IIIIl Specifically, if information is com- 
bined from xmi « 0.06 - 0.70 (7-band, total 42 MHz 
bandwidth), the constraint can be significantly tightened, 
i.e. AJnl = 0.16 for Omniscope (10 times smaller than 
an ideal CMB experiment), and Af^L — 2.8 by SKA 
(two times smaller than Planck). A prior of Ares = 
0.014 from Planck+WMAP CMB measurements [loj| cor- 
responds to a prior of A^esmr ~ 39, much larger than 
allowed by SKA and Omniscope alone, so adding this Tes 
prior cannot improve A/jvl from these experiments. Re- 
sults above assume linear perturbation theory is accurate 
to /cmax = 0.15 Mpc~^. If wc take a more conservative 
upper limit, fcmax = O.lOMpc"^, then AJnl of Omnis- 
cope is relaxed by a factor of ~ 2, i.e. AJml = 0.36, still 
4.4 times smaller than an ideal CMB experiment. 

Comparison with previous work. — Previously, Ref. [l^ 
reported forecasts of Af^L = [100,700,50,4] for 
[MWA512, LOFAR, SKA, Omniscope] based on informa- 
tion from a single redshift bin at the 50%-ionized epoch. 
Their results can be compared to the bottom set in Ta- 
ble llll Our results differ for a number of reasons: (1) 
Ref. [ij] computed the scale-dependent signature of PNG 
in Pat by fitting their approximate semi-numerical sim- 
ulations of reionization. In Ref.[3|], we showed by our 
analytical derivation and numerical LPTR calculations 
that this fit underestimates the scale-dependent bias due 
to PNG significantly. These differences are reflected in 
the smaller AJml we obtain for LOFAR, SKA, and Om- 
niscope. (2) The anticipated MWA512 configuration as- 
sumed by Ref. IJ] was also overly optimistic, while we 
adopt the current MWA128 configuration [s^. 

To model the effect of PNG on the 21cm power spec- 
trum, Ref. [ll] assumed the simple functional dependence 
of ionized fraction on local overdensity, used for illustra- 
tive purposes by Ref. [42| for the Gaussian case, to derive 
an ionized fraction bias bx for PNG. Unfortunately, they 
incorrectly used b^ to relate the 21cm power spectrum to 
the matter power spectrum (see Ref. [3|). 
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Conclusions. — This Letter suggests two approaches to 
constrain f^^ with 21cni power spectra from the EOR. 
If we take a conservative approach, i.e. assuming nothing 
about xiii{z) and b'^^^^{z), then the Pheno-model can be 
employed for a single redshift bin to provide Jnl con- 
straints with the same precision as a reionization model 
in which the reionization history is uniquely specified by 
a set of model parameters. However, using the ESMR, 
we demonstrate that a well-motivated reionization model 
can improve A/at^ in two ways: (1) A pathfinder mea- 
surement at a single redshift can best-fit the values of 
model parameters, which then can be used to estimate 
the desired redshift corresponding to xhii ~ 0.25 — 0.3, 
i.e. the "sweet spot" for cosmic-variance-dominated ex- 
periments. This can help tune single-epoch observations 
for maximum precision. (2) Multi-epoch measurements 
can be combined to improve A/jvl. We find that multi- 
frequency observation by SKA can achieve A/nl ~ 3, 
providing a new method to constrain PNG independent 
of CMB measurements, but with a precision comparable 
to PlancKs. A cosmic-variance-limited array of this size 
like Omniscope can achieve A/jv^ ^0.2, improving cur- 
rent constraints by an order of magnitude. These high 
precision observations may someday shed light on infla- 
tionary models. 
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